Weak Coulomb blockade effect in quantum dots 
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We develop the general non-equilibrium theory of transport through a quantum dot, including 
Coulomb Blockade effects via a 1/N expansion, where N is the number of scattering channels. At 
lowest order we recover the Landauer formula for the current plus a self-consistent equation for the 
dot potential. We obtain the leading corrections and compare with earlier approaches. Finally, 
we show that to leading and next leading order in 1/N there is no interaction correction to the 
weak localization, in contrast to previous theories, but consistent with experiments by Huibers et 
al. [Phys. Rev. Lett. 81, 1917 (1998)], where N = 4. 
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The Landauer formula relates the two-terminal con- 
ductance of a mesoscopic conductor to its quantum me- 
chanical scattering matrix [J]. This formula, and its ex- 
tensions to multiterminal geometries and other transport 
properties, has become one of the cornerstones of the the- 
ory of quantum transport. However, because the Lan- 
dauer formula relies on the single-particle scattering ma- 
trix, effects of electron-electron interaction are neglected. 

In a chaotic quantum dot the dominant part of the 
interaction is 0,13 

H c = E c (N d -N) 2 , (1) 

where E c = e 2 /2C is the one-electron charging energy in 
terms of the dot's capacitance C, is the number of 
electrons on the quantum dot, and J\f is an offset deter- 
mined by the voltages on nearby gates. The simplest way 
to include this interaction into the scattering formalism 
is to replace H c by a self-consistent potential 



(2) 



where the expectation value of the number of electrons 
{Nd) is computed in the presence of the potential Va- The 
self-consistent potential is of particular interest for non- 
linear and time-dependent transport Q. It is generally 
believed that Eq. is valid when the charge on the dot 
fluctuates freely either due to large contact conductance 
or high temperatures, in which cases Coulomb blockade 
is suppressed |2|. However, in this limit no systematic 
evaluation of the corrections to the quantum transport 
has been given for general non-equilibrium conditions. 
In this Letter, we develop a theory of non-equilibrium 
transport that reproduces the Landauer formula with the 
self-consistent potential at lowest order and admits cor- 
rections as a power series in the inverse of the number of 
scattering channels N. 

The situation we study is thus referred to as 'weak 
Coulomb blockade' 0], as a hard Coulomb gap in the 
current- voltage characteristic is not present. Histori- 
cally, the influence of quantum charge fluctuations on the 
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FIG. 1: Left: Schematic drawing of a quantum dot coupled 
to source and drain leads via point contacts. Right, top: a 
fictitious lead is side-coupled to electrons leaving the dot as 
a mathematical trick. The scattering Hamiltonian Hi, when 
included to all orders, decouples the fictitious leads, as illus- 
trated in the bottom right panel. 



transport properties of mesoscopic junctions first focused 
on the coupling to the 'electromagnetic environment' of 
an electron moving through a single junction This 
situation was recently extended to a single open contact 
in large-iV limit, similar in spirit to the present paper, 
but restricted to a point scatterer Q. To deal with the 
strong quantum charge fluctuations of almost fully trans- 
mitting point contacts, Flensberg and Matveev_ 
introduced a bosonized description (see also Ref. 
The weak Coulomb blockade situation has also been ad- 
dressed starting from the tunnel junction limit, using the 
phase as variable 0, and good agreement with exper- 
iments on large conductance metallic systems has been 
achieved Recently, the technique was extended to 

junctions of arbitrary transparency |12j . A crucial simpli- 
fication of all these approaches was that once an electron 
enters the dot, it never coherently returns to the leads. 
This was assumed to be satisfactory for large dots where 
the typical dwell time ra by far exceeds the time scale set 
by inelastic processes. 

For the problem we consider here, the coherent return 
of the electrons to the point contacts is crucial. Only 
then is transport in the absence of interactions described 
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by a unitary scattering matrix S. For the coherent case, 
the interaction correction to the dc conductance of an 
open quantum dot was first calculated in the bosoniza- 
tion [3, Il3l and later in the phase representation for- 
malisms |l4|. The results of these two approaches are 
strikingly different. Our secondary aim in this Letter is 
to clarify this situation. In this regard, we find that the 
result of Refs. 0, El > which is a formal expansion in the 
scattering matrix S, is not applicable for unitary S, de- 
spite claims to the contrary. However, it is applicable 
for a subunitary scattering matrix, which may serve as 
a phenomenological description of a quantum dot with 
strong relaxation. On the other hand, the 1/N correc- 
tions we obtain here are fully consistent with the analysis 
of Ref. 0| . Furthermore, whereas Ref. intended to 
describe the case of metallic quantum dots, i.e., neglect- 
ing quantum corrections to the non-interacting conduc- 
tance, we find that the same expression for the conduc- 
tance also describes the interaction correction to quan- 
tum interference corrections. Thus our theory applies to 
experiments on the weak localization corrections to con- 
ductance |l5j |. 

We now turn to a description of our calculation. We 
consider a quantum dot coupled to source and drain leads 
through point contacts with Ni and N 2 channels respec- 
tively (iVi + N 2 = N), see Fig.[T] The starting point of 
our calculation is the one-dimensional description of the 
electrons in the point contacts 

N 



3=1 



dx 



(,:) 



Here vf is the Fermi velocity (differing velocities in each 
channel may be removed by rescaling), rp^j(x) and tpuj 
describe left and right moving fermions in channel j, H c is 
given by Eq. (JJJ , and we set h = 1 . The lead-dot interface 
is taken at x = 0, the region x < corresponding to the 
leads. Following Ref. [3, we write the number of electrons 
in the dot as 



N d = N - dx: W Lj (x)^ Lj (x)+^ Rj (xW Rj (x)} : (4) 
J —00 

where No, the combined number of electrons in the dot 
and the leads, is taken to be a constant. 

The fields ip L and ip R in Eq. are not independent, 
since they are connected by the dot. After the substitu- 
tion which removes interactions inside the dot, their 
relationship is described by the scattering matrix 

$ Li (V,t)=Yl I dt ' Sy(M')W7'0, (5) 



where rj is a positive infinitesimal. (Scattering from a 
quantum dot with time-dependent potentials is described 
by a scattering matrix with two time arguments, see e.g., 
Ref. H3.) 



Instead of solving the interacting Hamiltonian J3J with 
the boundary condition J3J), we reformulate the problem 
using a trick that allows a treatment with standard per- 
turbative techni ques . Our method is similar in spirit to 
that of Refs. pi. Il3l Il6j |. As shown in Fig. 1, we cou- 
ple a fictitious reservoir to the left-moving electrons in 
the point contact. Then the problem reduces to the ex- 
actly solvable theory of Flensberg and Matveev, where 
electrons that enter the quantum dot never return coher- 
ently. We restore the original model by adding a term 



H x = 



N 

2v F J2 ($y (0)^ (V) + $ Li (v)i> Lj (0)) , (6) 



3=1 



to the one-dimensional Hamiltonian Eq. thus shunt- 
ing the fictitious leads and regaining coherent scattering 
by the quantum dot. 

We proceed by evaluating the current in perturbation 
theory in H\ using the Keldysh formalism and Eq. JSJl, 



I(t) = (TJ(t)e^^ dt '" l( - t ' ) 



(7) 



Here T c denotes ordering along the Keldysh contour c 
and the average is taken with respect to the Hamiltonian 
Ho . The current I is calculated as the weighted difference 
of the currents in the two point contacts, 

I = -ev F J2^n Ub^R^~^l^lM > (8) 



where A, 



Sij{N 2 /N) for j = l,...,JVi and A, 



(3) -SijiNt/N) for j = N 1 + l, 



,N. 



For the perturbative evaluation of Eq. J7J) we require 
the correlation function lal 



Tc n V s ! 4p (*p)^l i v k (4) 

f(t p ,t' q )f(s p ,sX 1 A 



p.q 



n 



f(tp, S' q )f(s p , t' q ) 



(9) 



(here all fields are taken at x — 0) where (• • • ) n j repre- 
sents the corresponding correlation function of the non- 
interacting theory (with E c — ► 0) and 



ln/(t, t') = const. — 



E n N 



dC,e 



-E c N^/tt 



(10) 



x In 



sinh [irT {t - t' + C - zO+sgn c (t))] 
sinh [ttTC] ' 



Here sgn c (i) is ±1 for the forward and backward branches 
of the Keldysh contour respectively. Additional 'anoma- 
lous' correlators containing N ip L and N ip R (or vice 
versa) give oscillating contributions to the current. These 
appear at Af th order in perturbation theory, so they are 
not captured by 1/N methods 3]. 
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Expanding Eq. J7J| to second order in Hi is now a straightforward calculation. The result is 



I(t) 



-trA/iR 



icvf 



dt'ds'tv [AS{t, t')G%(t', s')S f {s', *)] 



2tt ^ 2 

-evplm J dsdt'ds'tr [AG%(s,t)S(t,t')G%(t',s')S^(s',s)] sin[K (M')] 

I 



f(t,f)f(s,s') 



/(*, B')fM 



1/N 



(11) 



<t>: 




FIG. 2: a) At leading order in 1/N, we retain only tree 
diagrams, b) The one-loop diagrams 



where we wrote NKo(t, t')sgn c (t) = Imln /(t, t'), 



i.e. 



(t, = ^(l-e-^'-'J/j^'-f). (12) 



In Eq. (| 1 1|) the trace is over the scattering channels j 

r Rij { 



1, . . . , AT. The Keldysh Green's function GjL(i, a) is 



J 1 g— t f* dafirti(o-) 
GRij(t, S) — 5ijP . — , 

J vf sinn [ttI (t — s)\ 



(13) 



where "P" denotes the Cauchy principal value and is 
the (electro-)chemical potential of electrons entering the 
dot through through the ith channel, i = 1, . . . , N. With 
zero chemical potential for the left movers, electroneu- 
trality demands Y] j fiRj — 0. A similar expression holds 
for G\, but without the phase factor in the exponent. 

The linear dc conductance of Refs. is recov- 

ered upon taking fim — —eV s d^ii> i = 1,. ■ ■ ,N, setting 
S(t,t') = S(t — t'), and expanding Eq. Ijllfl in the source- 
drain voltage V s d- Truncating the perturbation theory 
in Hi at second order corresponds to a formal expan- 
sion up to second order in the scattering matrix S - 
which is precisely what was done in Refs. 

HE How- 
ever, as 5 is unitary, we find no justification for stopping 
here. That this approach gives unphysical results may be 
readily seen by substituting into Eq. the scattering 
matrix of a ballistic quantum wire with no backscatter- 
ing, resulting in a spurious correction to the Landauer 
conductance of e 2 /h per spin channel. 



In order to do a calculation to all orders in Hi, we 
need a different principle to organize the perturbation 
theory. In fact, the correlation function © shows that 
a systematic expansion in 1/N is possible [keeping EqN 
in Eq. (|12|l constant], while including Hi to all orders. 
The result is a sum over partitions into 'bubbles' in stan- 
dard diagrammatic language. Each bubble is formally of 
order N, as it involves a trace over the channel index. 
In order to be non-vanishing, bubbles must be connected 
by 'interaction lines' A^" 1 In f(t, t'). Denoting these by 
dashed arrows, as in Fig. the arrowhead to indicate 
the unprimed time, we see that at leading order we must 
have precisely one factor A -1 In f(t, t') for every bubble, 
except for the bubble that contains the current operator. 
The result is a tree structure, as shown in Fig. We 
thus find 



I(t) 



e 

ievF 



tvA^R (14) 
dt'ds' tiAS(t, f)G%(t', s')S^(s', t). 



Here G% (t',s') = Gf (*', S ')e-**V) where 4>(t'i,f 2 ) sat- 
isfies the self-consistent equation (see Fig. Et), 



<j){t\ s') = 2 / dtj c (t) [Ko(t, t') - Ko(t, s')} , (15) 
Gf(t,t) 



3c(t) = -^tr 



dt[ds[S(t, t'i)G% (t'i, «i)5t(«i, t) 



(16) 



The current j c is the particle current leaving the dot 
due to scattering inside. The quantity <fi(t', s') may 
be thought of as the difference in phases accumulated 
by electrons that enter the dot at times t' and s' in 

the presence of a Coulomb potential. Indeed, writing 

t' - 
cj)(t',s') = J , dteVd(t) and using Eq. fT2"|) for k , the 

self-consistency relation l|15|l takes the simple form 



Cd t V d (t) = ej c {t) ~ {Ne 2 /2n)V d (t), 



(17) 



which is precisely the equation for the dot potential in 
the self-consistent treatment of the interaction Hamil- 
tonian (JTJ. Hence, evaluating the current at tree level, 
we recover precisely the Landauer formula Eq. (|14J) . in- 
cluding the approximation (J2J for the interaction 0]. 
The leading correction to the above self-consistent the- 
ory adds one more 'interaction line' to each tree, which 
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results in diagrams with one loop, and generates the 
RPA-like series in Fig. [5Jd. The summation of these 



leads to a renormalization of the interaction line Eq. i|12|l : 

«(*,*') = «o(*,t') + f dtix{t,h)K(h,l/), where 



The four one-loop diagrams are shown in Fig. [2Jd, where the solid interaction lines denote n(t,t') [19) . We thus find 
the leading interaction correction to the current is 

7i(t) = ev|lm /" dti<dsid^d4K(ti,4)trS , (ti,ti)Gf (i' 1 ,si)S ,t (s / 1 ,f)AS , (f,4)G , f (4, 4)^(4, fi), 

leading to our final result for the one-loop correction to the dc conductance 

ott~j J cosn" [e/zi'j 

r 



ST , T Im J dsd J &U ^ 2 + ^ T K (co)tr [AS(e)AS\e + u)-AS(e)S\e + co)S(e)AS\e)] . (18) 



The result (|T%|) coincides with that found in Ref. |14j . 
We note that it gives not only the leading interaction cor- 
rection to the semiclassical conductance, but also to the 
quantum interference corrections. In particular, an en- 
semble average of Eq. I|18[) predicts the absence of an in- 
teraction correction to weak localization, which is in dis- 
agreement with the conclusions of Ref. [13, but consistent 
with experiments by Huibers et ai, who find that weak 
localization is well described by the non-interacting the- 
ory, suitably modified for the presence of a small amount 
of dephasing, down to the lowest temperatures available 
[Isj ]. There is a simple reason for this. The interaction 
correction Eq. Ijl8(l derives from a renormalized scatter- 
ing matrix S'(s) = S(e) + 5S(e), with 0O3 

6S(e) = -if f^ K M I tanh K e - w ) / 27 1 S ( £ - w ) 

+ tanh [(e + u) /2T] S(e)S*(e + u)S(e)] .(19) 

This correction preserves, at each energy, the circular 
invariance of the ensemble of scattering matrices that 
describes open dots. Quantities like the weak localiza- 
tion correction, that do not depend on correlations at 
different energies, are thus completely unaffected by the 
replacement S(e) — > S'(e). Note that this is a statment 
about the ensemble averaged conductance, and the con- 
ductance of any particualar dot configuration will ac- 
quire a temperature-dependent renormalization of ran- 
dom sign. We note that for non-ideal contacts with reflec- 
tion matrices, r C j, j — 1,2, Eq. (|T%)) predicts the renor- 
malization of the classical conductance g = 5152/(31+32) 

S 9j = -2^^,(1-^)1 In ^g^, rj 1 « T « gE c 
9j = — txr c ^r\- is the conductance of each contact 0. 



We have presented a formalism that elucidates the 
relationship between the conflicting theories of weak 
Coulomb blockade 0, 0, Q|. Our formalism provides 
systematic corrections to the use of a self-consistent po- 
tential in the scattering approach to quantum transport. 
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by the Packard Foundation and by the Danish Natural 
Science Reseach Council. 
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